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Abstract. In this work, we present a proof of the existence of real and ordered 
solutions to the generalized Bethe Ansatz equations for the one dimensional Hub- 
bard model on a finite lattice, with periodic boundary conditions. The existence 
of a continuous set of solutions extending from any U > to U = cxd is also shown. 
This continuity property, when combined with the proof that the wavcfunction 
obtained with the generalized Bethe Ansatz is normalizable, is relevant to the 
question of whether or not the solution gives us the ground state of the finite 
system, as suggested by Lieb and Wu. Lastly, for the absolute ground state at 
half-filling, we show that the solution converges to a distribution in the thermo- 
dynamic limit. This limit distribution satisfies the integral equations that led to 
the well known solution of the ID Hubbard model. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.30.+h. 

1. Introduction 

In 1968, Lieb and Wu I LWT] solved the one-dimensional Hubbard model H , and showed 
the absence of Mott transition in the ground state. A new version of this result |LW2j . 
with additional details on the exact solution for the ground state, was recently published, 
motivated by the increasing importance of the ID Hubbard model in Condensed Matter 
Physics. 

Despite being one of the most studied models to describe interacting electron systems, 
many questions concerning the Hubbard model remain as open problems, and the ID 
Hubbard model is the only one for which an exact solution was found. Still, many results 
concerning the Bethe Ansatz to the Hubbard model rely on unproved hypotheses. The 
purpose of this paper is to answer some of the questions raised in the recent paper by Lieb 
and Wu, concerning the existence of a solution to the Bethe Ansatz equations for the ID 
Hubbard model on a finite lattice, and whether or not this solution gives us the ground 
state. We will also discuss the solution in the thermodynamic limit. 

Consider a system with N electrons on a lattice with Na sites. The one-dimensional 
Hubbard model is described by the hamiltonian 

H = - ^^{c\^l^„c^^„ + h.c.) +U^ni^^nil, 
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where cj^, Cj^o- are the fermion creation and annihilation operators for an electron with 
spin a =t or J,, at site and nj-f, are the corresponding occupation number operators. 
We are using here periodic boundary conditions {Na + 1 = 1), so we are actually working 
on a ring with A^a sites. Particle- hole symmetry permits us to restrict to the case iV ^ Na- 

Since the hamiltonian commutes with X^i^i ^'^'^ 'Y^=\'^ih look for energy 

eigenstate with a fixed number of spins up and down. Let us denote these states by 
|M, M' >, where M is the number of spin-down, and M' the number of spin- up electrons. 
Spin-up and spin-down symmetry allows us to restrict to the case M ^ M' . 

We can expand \M,M' > in a basis of states of localized electrons, 

\M,M'>= f{xi,...,XN)\xi,...,XN >, (1.1) 

1 ^ a-, ^ Na 

where \xi, . . . , xj\f > denotes the state where the spin-down electrons occupy the sites 
xi, . . . , xm, and the spin-up electrons are located at xm+i, ■ ■ ■ ,xn- 

Our states are defined in the region R = {1 ^ Xi ^ Na , i = 1, . . . , N}. Given any 
permutation Q : {1, . . . , N} — > {Ql, . . . , QN}, we can define the region 

Rq = {1 ^ XQi ^ XQ2 ^ ... ^ XQN ^ Na} C R 

It is clear that 

UqRq = R, 

where the union is taken over the A'^! possible permutations. We should always keep in 
mind that the xi, . . . ,xm denote the positions of the spin-down electrons. 

The solution of the ID Hubbard model is based on a generalized version of the Bethe- 
Ansatz |Bej . For a sequence of unequal and ordered numbers ki < k2 < ■ • • < kj\f, we 
assume that the energy eigenstates are given by where the amplitudes f{xi, . . . ,xj\f) 

are given in the region Rq by 

N 

f{xu...,XN) = Y.AQ{P}e^p[tY,kp,Xj]. (1.2) 
p j=i 

The (A^!)^ coefficients Aq{P} have to be determined in order for the state to be well 
defined in the boundary of two regions Rq and Rq'. Also, we should make sure that the 
resulting state is antisymmetric with respect to exchange of two identical particles, and 
that it satisfies the periodic boundary conditions. 

Those restrictions result in a system of (A^!)^ linear and homogeneous equations for the 
coefficients Aq{P}. Fortunately, these equation were studied in detail by Gaudin [U| and 
Yang 0, in the context of a system of fermions with delta interaction. The only difference 
here is that we are working on a lattice, and not in the continuum, so the kj that appear 
on the system of equations have to be substituted by sin kj, which does not interfere with 
the algebraic analysis. 

The condition obtained in and p| for the existence of a non-trivial solution for the 
coefficients Aq{P} can be written as the generalized Bethe-Ansatz equations 

I sm kj — lAp + (7/4 

1^ isinkj-iAa-U/A ^ _ 1^ -iAp + iA« + U/2 ^ m n A\ 

j=l ■' /3=1 
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in terms of our original set {kj}jLi, and a set of auxiliary parameters {Aa}aLi- 
Defining 

e{x) = -2ta.n-^{2x/U) 
and taking the logarithm of the equations above, we obtain the Lieb-Wu equations 

Nakj = 2nl, + Y ^(2 sin kj - 2Ap), j = l,...,N , (1.5) 

/3 

^^(2sinA;, -2A,) = 27rJ«-^^(A«-A^), a = l,...,M. (1.6) 
i /3 
The coefficients Ij are integers if M is even, and half-integers otherwise. Similarly, Jq, are 
integers if — M = M' is odd and half-integers otherwise. We will restrict to the case 
where M is odd, and A^ is even, since in this case the ground state is unique for every U. 
For the ground state, we can choose the Ij and Jq, that will give us the correct solution in 
the limit U — > oo. In this limit, the equations decouple, and we have 

- ■ ^ + 1 _ M + 1 

^3 — 3 2 , Ja — Oi 2 • 

In the thermodynamic limit, taking A^q, A^, M, M' oo, keeping their ratios fixed, H1.5|) 
and ()1.6() become integral equations. The analysis of the resulting equations led to the 
well known solution of the ID Hubbard model. 

Some questions concerning the finite lattice case remained as open problems. The 
purpose of this paper is to provide rigorous results to some of those problems. In particular, 
we are going to show that the equations H1.5() . H1.6() do indeed have a solution, and that 
it is ordered in j and a (fcj+i > kj, Aa+i > A^). The ordering is important in the 
derivation of the integral equations in the thermodynamic limit. We also show that there 
is a continuous curve of solutions extending from [/ = oo to any U > 0. We are ready to 
state our main result: 

Theorem 1.1. The equations ()1.5() and (|1.6() have at least one real solution satisfying 

— TT ^ ki < ■ ■ ■ < kj\[ ^ -K, 

Ai < • • • < Am. 

Furthermore, there is a continuous curve of solutions defined by 7(t) : [0, 1] — > M''^^*^'''-^, 
where 

7(t) = {ki{t),...,kn{t),Ai{t),...,AM{t),U{t)), 

extending from any U > to U = oo . 

We know that the state obtained by the solution to H1.5() . ()1.6|) correspond to the true 
ground state in the limit U = oo, since we chose the coefficients Ij and Jq, in order to 
match these states. To prove that the ground state is given by the Bethe Ansatz solution 
for any positive U, we need a few additional facts. Since M and M' are odd, the ground 
state is non-degenerate for any U |LW2j . and our result shows that starting with the 
solution at f/ = oo, and going along the curve 7(t), the state is indeed the ground state, 
provided f{xi, . . . ,xn) defined by (|1.2|) does not vanish for all {xi, . . . ,xn)- 

On section 4, we prove that the generalized Bethe Ansatz equations have a solution 
that is algebraic in U . Since the Bethe Ansatz wavefunction is a rational function in the 
variables Zj = e^^^ , Aq and U (see |Wj). its norm will have at most finitely many zeros 
and poles as a function of U . We show that we can redefine our wavefunction at these 
problematic points, in order to obtain the physical solution for all U . 
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Once the existence of a solution that corresponds to the ground state of our finite 
system is shown, our next goal is to determine what happens with the set of k-s and A-s 
in the thermodynamic limit. As we increase the lattice size, a sequence of distributions 
measuring the density of A;-s and A-s can be defined. In the thermodynamic limit, we show 
that a subsequence of these distributions converge to the solution of the integral equations 
derived by Lieb and Wu. 

We are particularly interested in the absolute ground state at half-filling, in which case 
an exact solution can be found, including an explicit expression for the ground state energy. 
In this case, we show that the whole sequence converges to the explicit solution obtained 
in pTT] . 

2. Existence of a solution 

The goal in this section is to prove the existence of solutions to the Lieb-Wu equations 
that are ordered in j and a {kj^i > kj, A^^i > A^). The proof will follow from Brouwer's 
fixed point theorem, but first we need to define a convenient map to apply the theorem. 

2.1. Defining tiie map. Because of the symmetry of the coefficients Ij and Ja, we can 
look for symmetric solutions of the form kj = — /cat-j+i, and Aq, = — Ajv/-o+i- Therefore, 
our set of k-s and A-s are totally defined by 

N 

k = {kj} , j = - + !,..., N , 

TA 1 M + 3 

A = {A„} , a = ^— ,...,M . 
The other components of k and A, are obtained by the symmetry condition. In particular, 

^(M+l)/2 = 0. 

N , M-1 JV I A/-1 

We can define the map (p ■.M.2~'~ 2 ^1^2^^ 2 by 

<A(k,A) = (k',A'). 

where A' is defined by 

^9{2smk, - 2A'J = 27rJ„ - ^^^(A, - A^), (2.1) 
j P 
for {M + 3)/2 ^ a ^ M, and k' by 

k'^,^ =max{0,-^(2^/iv + V0(2sin/c^ -2A^))}, (2.2) 

a p 

whereas for i > ^ + 1 

k', = max{A:,Li, ^(2^/. + ^0(2 sin A:^ - 2A^))}. (2.3) 

" /3 

The function A'(k, A) is well defined since the left side of (|2.H) is continuous and strictly 
increasing in A^. It is also clear that k'(k. A) is continuous, so the map (/)(k. A) = (k', A') 
is continuous. 

It is also important to observe that 



A'q, < Cn,u, 
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where Cm,u is a constant for each value of N and U . Indeed, from (|2.1j) we have 
^6'(2sinA;j - 2A'„) ^ ~ 27r + Mtt ^ {N-1)tt, 



which imphes 



4A' 



a 



TT vr 



< tan - - — + 



U V2 2NJ [/' 

or 

U /"IT TT 



A'„< jtan(--— = 



2.2. The fixed point theorem. Now we are ready to state the result that allows us to 
use Brouwer's fixed point theorem. 



AT I M-1 



Lemma 2.1. If we define the domain 17 C M 2 + 2 hy 

^ ^Jv , 1 ^ ... ^ ^ TT, 

2 

^ Aiu+3 ^ ... ^ Am ^ Cn,u, 
2 

we have 

(p{n) c n. 

Proof. By definition, 



^ ^tLj^t, ^ • • • ^ ^N- 



Also, 



. N -I 

2tt 



N 

< TT^-- < vr. 



Regarding A^, we have 

6(2 sin kj - 2A'J > ^ A'„ > 0. 

Furthermore, the right side of (|2.1j) is increasing in Aq and Ja+i > -/a, whereas the left 
side is increasing in A^. Restricting the map <j) to this implies 

Combining these results with A^j < Cn,u we conclude the proof of the lemma. 

Since O is compact and convex (in particular, it has no "holes"), there is a continuous 
and invertible map 

where D^^ ~ is a disk in (A^ + M — l)/2 dimensions. The map 

_1 iV I M-1 N I A/-1 

IpOCpOlp : L) 2 2 -^D2~^ 2 

is continuous and takes the disk into the disk, and by Brouwer's fixed point theorem, it 
must have a fixed point. Since "0 is one-to-one, the fixed point of '0 corresponds to a fixed 
point of (p. 
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2.3. The fixed point is a solution to the Bethe Ansatz equations. Since we ar- 
tificially introduced the maximum in the definiton of k', we need to show that the fixed 
point (k*, A*), such that 

(/.(k*,A*) = (k*,A*), 

is indeed a solution to the original equations. 
In other words, we need to show that 



<k*N < ■■■ <k 



Taking kj = 0, we have 



2-kI- 

k'. = > ^ k\^^ > 0. 

■' jyJa 2+^ 



That also implies 



k*N+i = -TT [2vr/iv^i + ^ e{2 sin k*^^^ - 2A^) 



1 

Let us consider now 



= max{/c^^^, -^[27r/^^3 + ^ 0(2 sin A:^_^3 - 2A^)] }. 

2 2 J V /-I 2 2 

a ^ 

Assuming that A;^ „ = k*f^ -, , we have 

-2-+-i "2"+l 

^jv , n = max-! kff ,ki^! -\ — — — [ > /cat -. , 
which is a contradiction. Therefore, k% > k*^^ , and proceeding by induction we get 

-2-+3 "2" + -'- 

kN__^^ < ^iV_|_3 < • • • < kj^. 

We have therefore concluded the proof that the fixed point (k*,A*) satisfies H1.5() and 

dm). 

3. Existence of a continuous curve of solutions 

We want to analyze now what happens to our solutions as U changes. Our goal is to 
show that there is a connected set of solutions extending from any U > to U = oo. 

From what we have seen in the last section, we know that on the part of the boundary 
of fl defined by kN,, = or ki+i = ki, the vector field defined by 

2 

Vt;(k,A)=0^(k,A)-(k,A) 

has a normal component that always points inward. Also, for = we have A^ > 0, 
and Aa = Aq,+i implies A^ < A.'^^^. Therefore, the vector field actually points inward in 
all the boundary of O. 

All we need to prove then is that given any family of continuous vector fields 

V, :D«^M" , iG[0,l], 

such that the normal component of points inward on the boundary of the disk D", 
there is a connected subset of x [0, 1] in which V((k, A) = 0, extending all the way 
from t = to t = 1. 

The proof of this statement is trivial if for each t € [0, 1] the vector field v is non- 
degenerate. (Here, we consider a vector field to be non-degenerate if, for every point x 
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such that v(x) = 0, the determinant of the jacobian \djVi\ does not vanish.) In this case, 
since the zeros of the vector field are always isolated, each zero at t = will follow a 
smooth curve as we increase t from to 1, by the implicit funtion theorem. 

In the general case, however, zeros can collide with each other and disappear, as we 
change t, so the proof is not that simple. Let us assume now that we have a general family 
of vector fields, possibly degenerate. 

We can always approximate by a polynomial vector field, in the sense that for every 
e > there is a polynomial vector field p such that ||vt — p|| < e. If is a polynomial, 
it will have a finite number of sets of connected zeros. We can define an index for every 
isolated set of zeros. For an isolated zero xq, we define its index |GPj by the degree of the 
map (j) : Se ^ 5""^ defined by 



where is a sphere center at xq with radius e, and e is sufficiently small so that xq is the 
only zero of in its interior. So the index is just a measure of how many times we map 
the sphere 5^ into the sphere S^^^. 

For any connected set of zeros, we can define the index analogously, by isolating the set 
from other zeros by a surface that can be smoothly deformed into a sphere. Alternatively, 
we could define the index by adding a small perturbation to remove the degeneracy, and 
adding the index of all resulting isolated zeros. A general definition of the index for vector 
fields can be found in |GSj . 

Since Vf always points inwards in the boundary of the disk, the sum of the indices of all 
connected sets of zeros is equal to one, for all t. As t varies from to 1, the zeros describe a 
continuous trajectory. A given zero of index +1 cannot simply disappear, unless it collides 
with a zero of index —1. Every connected set of zeros can at some point break up into 
other sets, provided the sum of indices is preserved. Also, new sets of zeros can be created, 
provided their total index is 0. 

Let Z C X [0, 1] be the set of all zeros of v, where we now drop the index t, and 
consider the vector field defined on the cylinder x [0, 1]. If a zero of vq with index 
+1 is not connected to a zero of vi by Z, then it has to be connected to another zero of 
Vq with index —I. Since the total index is always equal to 1, there are not enough zeros 
of index -1 to annihilate all the +1 zeros, and there will be at least one connected set of 
zeros from t = to t = 1. 

We can find a sequence of vector fields {w„}, such that v + w„ is a polynomial vector 
field and lim„_^oo |w„| = 0. For each n, we have a continuous curve of zeros of v + w„, 
from t = to t = 1. What we need to do is to prove that the curve cannot be disrupted 
as we take the limit |w„| — > 0. 

Let us also define the set 



where Dq is the disk at i = 0, and Di the disk at t = 1. Consider the topology r* in Z' 
defined by taking the intersection of the open sets in the metric topology t of D x [0, 1] 
with Z'. 

Let us assume for a moment that there is no subset of Z connecting Dq to Di. In this 
case, we can find two disjoint open sets Xq,Xi £ r* such that 



and Xq U Xi = Z' . But Xq and Xi are closed sets in the metric topology r (defined by 
the open balls in the usual metric on the cylinder), since they are given by zeros of a 




Z' = DqU ZU Di 



DqCXq 



DiCXi 



8 



PEDRO S. GOLDBAUM 



continuous vector field. Therefore we can find disjoint open sets Oi, O2 £ t such that 

XoCOo , XiCOi. 
But the set (Oi fl 02)*^ is nonempty and closed, so 

|v(x)|/o , Vx G (Oi n O2)', 

implies 

|v(x)|>(5 , VxE(0in02)'. 
Hence, for [w[ sufficiently small, 

[(v + w)(x)[ > 5/2, 

and there is no curve of zeros of v + w connecting Dq to Di . This is in contradiction with 
the existence of a sequence of polynomial vector fields converging to v. 

Therefore, we conclude the proof of the existence of a connected set of solutions of the 
Bethe-Ansatz equations, extending from any J7 > to arbitrarily large values. Along this 
curve, the energy given by 

N 

E{U{t)) = -2^ cos kj{t) 

is continuous in t. 

4. NON- VANISHING NORM OF THE WAVEFUNCTION 

As we discussed before, the existence of a continuous curve of solutions implies that 
the right side of is indeed the true ground state, provided that its norm is not zero. 
However, since we are working on a lattice, and not in the continuum, the wavefunction 
could in principle vanish (/(xi, . . . ,xn) = 0) , even though the coefficients Aq{P} 
are not all zero. Therefore, to complete the proof that (|1.H) is the true ground state, we 
need to prove that /(xi, . . . ,xn) does not vanish identically. 

We will follow here the same strategy used by C. N. Yang and C. P. Yang jYYj to 
prove that the Bethe Ansatz does indeed give the true ground state for the anisotropic 
Heisenberg model. 

Let us first define the variables Zj by 

f e^'^J , j = 1, . . . ,N, 
= \ Aj.N , j = N + l,...,N + M. 

We can write the equations H1.3|) . 1)1.4(1 as polynomial equations in the variables Zj and 
U. So we have a system of + M polynomial equations, in + M + 1 variables, 

Pi{zj,U)=0 , i = l,...,N + M, (4.1) 

and our goal is to find solutions given by algebraic functions Zj{U). 

Let us remind ourselves of some basic facts concerning polynomial equations in several 
variables. A more detailed background on the subject can be found in _CLO . Consider the 
space of polynomials in n complex variables, with complex coefficients, which we denote 
by C[xi, . . . , Xn]- Given a set of polynomials {/j} C C[xi, . . . , where i = 1, . . . , m, we 
can define an ideal /(/i, . . . , fm) C C[xi, . . . , Xn] by 

m 

I{fl,---,fm) = {^hifi , hi,...,hm& C[xi, . . . ,Xn]}- 
1=1 
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We can also define an affine variety V{I) C C" by 

= {(ai,...,a„) gC" : /(ai, . . . , a„) = , V/ G /}. 
V(I) is the set of solutions to the original polynomial equations 

/l(xi, ... ,Xn) = /2(X1, . . . ,X„) = • • • = fm{xi, . . . , X„) = 0. (4.2) 

We say that the set of polynomials {fj} form a basis to the ideal /. But an ideal can 
be defined by many bases. A particularly convenient basis for an ideal is the Groebner 
basis (see appendix A or |CLUj for details). The Groebner basis is a set of polynomials 
with the same roots as the original set, with the nice property that it contains polynomials 
where some of the variables are eliminated. We will be using a Groebner basis for the ideal 
defined by the Bethe Ansatz equations shortly, but we still need to define the projection 
of an affine variety. 

Given V{I), we can define vr^'^') : V{I) C C" ^ C by 

7r('')(ai, . . . ,a„) = Ofc 

This projection gives us the the values of Xk for which the system of polynomial equations 
has at least one solution. 

We are ready to state the technical result that will assist us in our proof. 

Lemma 4.1. Let V C C" be an affine variety. Then one of the following has to he true: 

• 7r(^)(y) consists of finitely many points (possibly zero), 

• C — -n^^^iy) consists of finitely many points (possibly zero). 

Proof: If we define the ideal /('^^ = /nC[xfe], given by the polynomials in / that depend 
only on the variable x^, we have 

where TV is a proper subset of V{I^^^). The right side of this relation is trivial, since 
(oi, . . . , Ofc, . . . , an) € V{I) implies that is a root of any polynomial in I^^). The left 
side follows from |GLOj (Theorem 3, p. 123). 

But y(/(*^)) and W are affine varieties in C, which can either be C or a set with finitely 
many points (roots of a polynomial equation p{xk) = 0). If V{I^^'^) = C, is a finite set 
and C — Tr^^\V) consists of finitely many points. Otherwise, tt^''\V) is finite. 

As a first application of this lemma, we see that our system of polynomial equations 
will actually have a solution for any complex U, except for at most finitely many values, 
since we proved that it has at least one solution for U real and positive,. If we choose the 
ordering zi > Z2 > ■ ■ ■ > zjy^M > U, our Groebner basis will not have any polynomial 
depending only on U. Therefore, our last polynomial in the basis should be of the form 

p{zM+N, U) = 0. 

There is nothing special about the ordering we chose, so can also construct a Groebner 
basis with the ordering zi > ■ ■ ■ > Zk > U , for any k between 1 and N + M. Therefore, 
the solutions of (|4.1|1 will have to satisfy 

p«(zi,C/) = 0, 
P^^\Z2,U) = 0, 



P^^'+'^HZM+N^U) = 0. 
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Let us assume for a moment that none of these polynomials vanishes identically. Then, 
we can factorize each of them into irreducible polynomials 

(k) 

Every solution will be given by roots of some combination of . Since there are finitely 
many combinations, and infinite solutions, there will be one particular combination that 
will give us solutions for all but finitely many values of U (we are again using Lemma 4.1). 

If one of these polynomials vanishes identically, it means that our system of equations 
is degenerate, and we can add equations to it, reducing the set of solutions, in such a way 
that the previous argument will work for the extended set of polynomial equations. 

Therefore, each of our variables Zk will be an algebraic function of U defined by 

pP^(zfc,C/) = 0. (4.3) 

Since we know that our wavefunction is a rational function in and U, if we go from 
J7 = oo to any finite U > along the set of solutions defined by ()4.3() . our wavefunction 
will have finitely many zeros and poles. Since those are isolated, we can redefine the 
wavefunction at these points by the limit of the normalized wavefunctions as we approach 
them, since the energy of those states varies continuously and is equal to the ground state 
energy for all U. Again, this idea goes back to |YYj . in the context of a spin system. That 
concludes the proof that the Bethe Ansatz gives us the true ground state of the system 
for any U > 0. 

If N/Na ^ 2/3, we could actually prove a slightly stronger version of this result, in 
which we show that the wavefunction that results from the solution of ()1.5|) and is 
not zero for all f/ > 0. 

Theorem 4.2. // the density of electrons per site satisfies 

N 2 

iv;; ^ 3' 

the norm of the state given by the Bethe-Ansatz method is strictly positive: 

II ^ f{xi,...,XN)\xi,...,XN > \\>0. 

1 ^X, ^ Na 

Therefore, the Bethe-Ansatz gives us the true and unique ground state of the system. 

Proof: Let us start by considering the case M = M' = N/2. Since we know that for a 
given solution of ()1.5() and H1.6() . the coefficients Aq{P} are not all equal to zero, we just 
need to show that /(xi, . . . ,xj\f) = in the region R implies 

Aq{P} = 0, 

for any P and Q. 

Let us consider the region Rq defined by 

Rq = {1 Xi ^ X2 . . . XN ^ Na} C R. 

Gaudin |U] showed that /(xi, . . . ,X]\f) in R is totally determined by its value in the 
region Rq, and Aq{P} = 0, for any P, implies Aq{P} = 0, for all P, Q. Therefore we just 
need to show that f{xi, . . . ,X]\f) = implies Aq{P} = 0, for any P. We will denote the 
restriction of /(xi, . . . , xn) to Rq by 

N 

fo{xi, ■ ■ ■ ,xn) Xj\ . 

p j=l 
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Also in [H], it is shown that for any two permutations P and Q, such that Q only acts 
on variables corresponding to one kind of spin, we have 

Ao{PQ} = I{Q)Ao{P}, 

where I{Q) is the sign of the permutation Q. This can be proved by combining the 
boundary conditions at the different regions Rq with the antisimmetry of the wavefunction 
with respect to exchange of identical particles. In particular, if P and P' differ only by 
the exchange of coordinates of two identical particles, 

Ao{P'} = -Ao{P} 

Therefore, we have 

/o(x2, xi, . . . , xn) = ^ Ao{P} exp [i{kpiX2 + kp2Xi H h kpNXN)] = 

p 

- ^^o{-P'}exp [i{kp/ixi + kp'2X2 H h kpNXw)] = -fo{xi,X2, . . .,xn)- 

P' 

Notice that this is not a trivial consequence of the antisymmetry alone, since /o is not 
the true wavefunction when X2 < xi. 

As a result, we have that /(xi, . . . ,xn) = implies fo{xi, . . . , xjy) = in the region 

{1 ^ Xi,...,XM ^ XM+1,---,XN ^ Na} C R 

If Na ^ 3N/2, that includes 

S = {1 ^ xi,...,XM < N,N + 1 ^ XM+i,---,XN ^ 3N/2}CR 
If /o is zero in 5, we can fix X2, • • • , xn and consider the system of equations 
fo{xi,X2, . . . ,xn) = , 1 ^ xi ^ iV, 

or 

e^*=i"iC«(x2,...,x,v) + --- + e'^^"iC7W(x2,...,x,v) = , 1 ^ xi ^ iV, (4.4) 
where each C^"'^ is given by a sum over permutations such that P{1) = n: 

= ^ AQ{P}QlCp[i{kp2X2 + ■■■ + kpNXN)\. 
P:Pl=n 

Therefore, 1)4.4(1 is a homogeneous system of N linear equations in the variables C^^\ 
The determinant 



does not vanish, since — vr < fci < • • • < fcjv < vr, and we have 

C^") = , 1 ^ n ^ iV. 
We can fix X3, . . . , xat, and write 

= ^ e*'^'"^2C("'™)(x3,...,X7v) = , 1 ^ X2 < 
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where 

Cin,m) ^ [i^kp^X^ + ■■■ + kpNXN)] • 

P:Pl=n,P2=m 

Again, the determinant does not vanish and 

After repeating this argument for all spin-down variables, we will have only AI = N/2 
pseudo-momenta /cj in each term. That is why it suffices to have iV+l ^ xm+It'-^xn ^ 
3N/2. By the time we get to the variable xat, we will have 

Ao{P} = 

for all P. 

If the set {ki} is given by a solution of the equations (|1.5j) and (|1.6() . they correspond 
to a non trivial set of coefficients Ao{P}, and by the result above, f{xi, . . . ,X]\f) 7^ for 
some {xi, . . . , X]y) £ R, and the state given by the right side of p.lj) has a strictly positive 
norm. That concludes the proof of the theorem. 

In general, M 7^ M', and the argument above still applies provided Na ^ N + M 
(which includes Na ^ 3N/2 particular case). 

5. The thermodynamic limit 

Let us now analyze the behavior of the solutions {/cj, A^} as we take the thermodynamic 
limit, by taking Na — > 00, keeping the ratios n = N/Na and m = M/Na fixed. 
We can define 

Pi^^^f = Wu. TT ' kj ^ k< kj+i, (5.1) 

and Pi{k) = for k < ki or k ^ k^- The above expression defines a sequence of functions 
Pi G L^([— vr, vr]). As i increases, so does the number of sites, and i — > 00 corresponds to 
Na 00. Analogously, we can define ai{A) G L^{R) by 

iVa I A/3+1 - 

and (Ti(A) = for A < Ai or A ^ Am- From the definition, the norm of pi and 
satisfy 

II II 1 II II M-1 1 
r 1 = = n , r 1 = = m . 

mwL , II *iiL 

In particular, these functions are uniformly bounded as we take the thermodynamic limit. 
Let us also define 

Mk) = Nak-Y,0{'^s^^k-2Ap). (5.3) 
/3 

From the Lieb-Wu equations, we have fi{kj) = 2711 j. Also, 

fm ^ Na + ^M ^ Na[l + ^^ 

which implies 

8m \ 

ir)' 

or 

{l + 8m/U) 



fiikj +1) - n{kj) = 27r ^ \kj+i - kj\Na{l + 
Pi{k) ^ 



27r 
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Since Pi{k) is bounded, pi £ LP([— vr, tt]), for 1 ^ p ^ oo. The sequence {pi} is actually 
uniformly bounded in L^. 

With a similar goal, we define 

5i(A) = ^(2 sin A:, - 2A) + 0(A - A^) , 

j 13 

such that gi{Aa) = 27r Ja. Its derivative is bounded by 

5^(A) ^ 77 , , 4.in.,-4 -A T2 ^ (5.4) 

i ^ + I V ) 

Once more, that implies 

8n 

gi{Aa+l) - gi{Aa) = 2tT ^ \Aa+l - Aa\Na — , 

or 

4n 

f^i(A) ^ -jj. 

nU 

The first inequality on (|5.4|) shows that all g'^ are bounded by a function in L^(M). There- 
fore, {(Tj} is uniformly bounded in L^. 

Since both {pi} and {cxj} are uniformly bounded, we can use Banach-Alaoglu theorem 
to show that there is a subsequence converging weakly in (see |LLj): 

Prii ^ Pw 5 ^ 

We will be working now with this subsequence, but we will keep the index as i just to 
simplify the notation. Our goal is to show that this subsequence also converges pointwise. 

If we take equation (|1.5() with index j, subtract it from the same equation with index 
j + 1, and divide the result by 2TrNa\kj-^-i — kj\, we obtain 

1 1 6'(2sinA;,+i - 2A«) - 6'(2sinA;,- - 2A/3) . , , , 

= S^-.^IAm-A,!. (5.5) 

where kj ^ k < kj^i and A^ ^ A < A^+i. We can fix k now and consider ()5.5() as an 
equation for pi{k). 

To show the existence of a pointwise limit, we just need to show that the right side of 
1)5. 5|) converges for all k. We know that the fraction on the right side converges to the 
derivative of 9, if \kj^i — kj\ — > 0. Since we did not prove that this happens for every k, let 
us divide the domain of p into Q € [— vr, vr], defined hy Q = {k € [— vr, vr] : j/cj+i — kj\ — > 0}, 
and its complimentary Q^. 

In Q^, the distance between the surrounding kj does not converge to zero. For a A; in 
the interior of Q^, we can find an e and a subsequence such that \kj — k\ > e/2, for all 
j. That implies that the weak limit vanishes in this interval of lenght e around k, since 
the integral of the characteristic function of any subset of this interval times pi converges 
to zero. Also, the subsequence will have a pointwise limit on this set that is equal to 
zero (almost everywhere). Since the boundary of Q'^ has measure zero, the value of the 
function at those points is not relevant in the context of L^-spaces, which are defined by 
equivalence classes of functions that differ only on a set of measure zero. 

Therefore, the subsequence {pi} has a pointwise limit in Q'^ that agrees with the weak 
limit. We will show that the same holds true in Q. A similar analysis can be done for 
the sequence {fjj}. We define the set i? G M by B = {A G M : |Aq,+i - A^I ^ 0}. In its 
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complimentary B^, the subsequence has a zero pointwise hmit that agrees with the weak 
hmit. 

Now if k G Q, the last term of (|5.5j) is 

'°'^.5^[0'(2sinfc - 2A^) +Ef A^)]a,(A)|A^+i - A^|, (5.6) 



13 



where £'|^'*(A;, A^) is the error in approximating the ratio in (|5.5|) by the derivative of 6. 
Since 9" is bounded, we have 

J]sf)(A:,A^)a,(A)|A^+i-A^| ^ ^ ~ = Cm\k,+, - k,\ ^ 0, 



/3 /3 



as z ^ OO. 

We need to show now that ()5.6() converges to an integral. We first notice that 
^6l'(2sinfe-2A^)cTi(A)|A^+i - A/3| = / 6l'(2sinfc - 2A/3)a,(A)dA + 

+ j [e[^{2s\iik- 2Kp) - 9'{2smk- 2A^)]o-i(A)(iA, (5.7) 

where 6[j. is the truncated version of 6', defined by a piecewise constant function, whose 
value on the interval [A^,A^+i) is equal to 0'(2sinA: — 2A^). Again, since the second 
derivative of 9 is bounded uniformly by a function in L^, we can write 

|C(2sinA;-2A^)-0'(2sinfc-2A^)| ^ g(A)[A^+i - A^|, 

where g{A) G L'^{R). It follows that 

lim / |6'i^(2sinA; - 2A^) - 6l'(2sinA; - 2A^)|cTi(A)(iA ^ lim / ^^dA = 0. 

Since {cjj} has a weak limit, the first integral on the right side of (|5.7)) converges, and 
we have 

lim pi{k) = — - [ 9' {2 sink- 2A)a^(A)dA, (5.8) 

j^oo 27r vr J £} 

for all € Q. A similar analysis for cTi(A) yields 

lim fji(A) = -- / 0'(2sin A: - 2A)p^(A)(iA + ^ [ 9' {A - A')a^{A)dA, (5.9) 
i^oo TT Jb 27r 

for ah A £ B. 

We just showed that the subsequences have a pointwise limit. But since {pi}, {cTj} are 
uniformly bounded, dominated convergence implied the existence of a strong limit, which 
equals the pointwise limit and the weak limit almost everywhere. Therefore, as we go to 
the thermodynamic limit on this subsequence, our k-s and A-s converge to distributions 
satisfying the integral equations 



p{k) = — + cosk [ K{sink - A)a{A)dA, 
2vr Jb 



(5.10) 



a{A) 



[ K{sink-A)p{k)dk- [ K^{A - A')a{A)dA, (5.11) 
JQ Jb 
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where we define 



K{A-A') = --^'(2A-2A') 



TT 



1 f°° 

K^{A-A') = 6'{A-A')= / K{A-x)K{x-A')dx. 

It was shown in |LW2j that it is convenient to restrict the integral in A; to a domain 
where the function sinfe is monotonic. Hence, we divide the range of integration Q into 
Q_ = Qn [-7r/2, 7r/2] and = Q n ([-vr, -tt/2) U {it/2, vr]). The first thing to show is 
that for every k = {it/2 + x) G Q+, there is a k = {it/2 — x) G 

Since the derivative of fi{k), as defined in (|5.3j) satisfies fl{k) ^ Na for k € [— 7r/2, 7r/2], 
the set Q- consists of a single interval [—k*,k*] (possibly [—tt/2, tt/2]). Since f{k) ^ 
f{k*) + Na{k — k*) for any k G [k* ,tt — k*], and the kj-s are ordered, we conclude that 
the existence of a gap in the interval [k* ,tt/2] implies that (5+ n [tt/2, tt — k*] = 0, and for 
every positive k in there is a corresponding point in Q_ given by the reflection with 
respect to A; = tt/2. The same holds for for every negative k in where now we need 
to take the reflection with respect to —tt/2. 

Therefore, if we define the set Qlj. by the reflection of Qj^ with respect to ±7r/2, 

Q; = {±(7r -k),keQ+r\ [tt/2, tt]}, 

and we will have C Q-. Finally, we can avoid the integration in Qj^ by noting that 
^^(sinA; — A) is an even function of A; — tt/2, which implies by (|5.1UI) that p{k) — \/2tt is 
odd as a function of A: — tt/2. Therefore, we can rewrite part of the first integral in (|5.11l) 
as 

/ + /" K{smk-A)p{k)dk= [ +[ ^ ~ ^) = A /" K(sinA:-A) 

Now sin A: is monotonic in the domain of the integral in k, and we can make the change 
in variables 

t{x) = ^ ^ , 

^ ' 2tt 
f{x) = (1 -a;2)-V2p(sin-ix), 

for —1 ^ X ^ 1. li A,D are the images of Q^,Q\ under sin A;, the integral equations 
(|5.1()|) . (|5.11|) can now be written 

/oo 
K{x - x')B{x')a{x')dx' , x e A, (5.12) 
-oo 

/oo 
K{x - x'){A{x')f{x') + D{x')t{x'))dx' 
-oo 

/oo 
K^{x-x)B{x)a{x')dx' , x e B, (5.13) 
-oo 



where A{x) = Xa\d^ B{x) = xb, and D{x) = 2xd {xx denotes the characteristic function 
of the set X). 

These equations have to be solved for x in the sets defined above. However, we can let 
the right side of (|5.12|) and (|5.13|) define f{x) and (j{x) for any x E M, by setting t{x) = 
for |x| > 1. Since the integrals only depend on f{x) and (j{x) at the original intervals, 
there is a one to one correpondence between the solutions to the extended equations and 
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the solutions of the original equations. In particular, the new equations will have a unique 
solution if and only if the solution to (|5.12j) and l|5.1,Sj) is unique. 
It is convenient to write the integral equations in the operator form 

/ = t + KBa, (5.14) 

a = KAf + km- k'^Ba, (5.15) 

where K is the convolution with K and A, B, D are the multiplication by A, B and D. 

The proof that (|5.14j) and (|5.15|) have a unique solution for given sets Q and B was 
done by Lieb and Wu in , LW2^ . We did not prove here that the sets Q and B are uniquely 
defined for every n = N/Na and m = M/Na, although we believe they are (and we also 
belive that they are given by intervals Q = [-kmax, kmax], B = [-Amax, -^-max])- However, 
this is not a problem for the absolute ground state of the half filled band, in which case 
an explicit solution for (|5.14)) and (|5.15|) can be found. In this case, we have n = 1 and 
m = 1/2, and Lieb and Wu proved that Q = [— vr, tt] and B = M. Let us restrict now to 
this case. 

Since the solution is unique, and every subsequence of {pi}, {cjj} contains a further 
subsequence that converges to the same p and a, we can actually prove that the whole 
sequence converges. Otherwise, we would be able to find a subsequence such that \\pn^ — 
pII > e or \\cTni — > e, which is a contradiction. 

The limit distribution in the thermodynamic limit for k-s and A-s are then given by the 
solution of (EHU, obtained in [iWT] . 

1 cos A; /"°° cos (wsin A;) Jo(a;) , 

^ ' 27: J Q cosh(a;C//4) ' ^ ' 

which yield the ground state energy 

E = -2Na r p{k) cos kdk = -ANa r t^'^^twl du;, 

where J„ is the Bessel function of order n. 

6. Conclusions 

We have concluded the proof that the Bethe Ansatz equations, and in particular the 
Lieb-Wu equations H1.5() and H1.6|l . have a solution for finite M, and Na- The solution 
is indeed real and ordered in the indices j and q, as assumed by Lieb and Wu in the 
derivation of the integral equations for the model in the thermodynamic limit. 

We have also shown the existence of a continuous set of solutions extending from any 
U > to [/ = oo. The continuity is important to show that the state obtained by the 
Lieb-Wu equations is indeed the ground state of the system, provided it does not vanish. 
Since we know that the state obtained is the genuine ground state of the system at U = oo, 
and that the ground state is nondegenerate for any U > (for M, M' odd), a normalizable 
state given by the Bethe Ansatz should also give us the ground state of the system for 
any U, since the energy cannot discontinuously jump from an excited level to the ground 
state as we increase the interaction strength from a finite value to the limit U = oo. 

To conclude the proof for a finite lattice, we show that the norm of the wavefunction 
obtained by the Bethe Ansatz equations is well defined and not zero, for all U > 0, except 
for finitely many values of U. If Uq is one of this values, the correct eigenfunction will be 
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given by the limit C/ — > C/q of the normahzed wavefunction, since Uq is an isolated zero or 
singularity. 

Going from the finite lattice to the thermodynamic limit, we focus primarily on the 
half- filled band N = Na- In this case, the set of k-s converge to a distribution p{k) on 
the interval [— vr, vr]. If we consider the absolute ground state with M = M' , the set of 
A-s converge to a distribution ct(A) in M, and we showed that p{k) and ct(A) satisfy the 
integral equations derived by Lieb and Wu, whose solution is given by (|5.16|) . (|5.17|) . For 
other fillings, we also prove the convergence (at least for a subsequence) of the distribution 
to the solution of the integral equations, but we do not prove uniqueness of the sets Q and 
B that contain the support of p and a. 

The author is indebted to Elliott H. Lieb for useful comments and suggestions. 

Appendix A. Groebner bases 

Every ideal /(/i, . . . , fm) generated by the polynomials {fk}^=i 7^ {0} has a Groebner 
basis. The general idea in constructing a Groebner basis is to eliminate one by one the 
variables of the polynomials fk, obtaining a set of polynomials G = {gk}k=i that generate 
the same ideal as {fk}k=i- One of its main applications is to solve systems of polynomial 
equations, since reducing the number of variables simplifies the task of finding solutions. 
This process can be thought of as a generalization of the process of solving a system of 
linear equations by elimination of variables. We will present here the essential facts about 
Groebner basis in order to understand the results of section 4. The interested reader 
should refer to |GLOj for a precise definition of a Groebner basis and the algorithm used 
to obtain it. 

Given an ideal / C C[xi, . . . , we can define the /c-elimination ideal by 

ik = In C[xk+i, . . . 

Ik is the set of all polynomials in / that do not depend on the first k variables Xk- What 
makes the Groebner basis special is that 

Gk = Gn C[xk+i, ...,Xn\ 

is also a basis for Ik- 

As a consequence of this property, if the variety defined by (|4.2j) consists of finitely 
many points, Gn-i 7^ {0}, and we can solve the corresponding equations for Xn, to find 
V{In-i)- Then, we can move on to Gn~2i and solve it for x„_2i finding V{In-2)i and so 
on. In the end, we will find V{I) just by solving polynomial equations in one variable. We 
are ordering the variables here by xi > X2 > • ■ ■ > Xn- This ordering determines which 
variables will be eliminated first. Of course we can choose any other ordering, and we will 
obtain a different Groebner basis, but the set of solutions should be the same 

As one example, consider the set of polynomial equations /i = + + z — 1 = 0, 
/2 = x + ?/ + 2; — 1 = and = x + y — z"^ — 1 = 0. Fixing the ordering x > y > z, the 
ideal /(/i, /2, /a) has G = {x + y + z — l,y'^ — y + yz — z, z"^ + z} as a possible Groebner 
basis. Therefore, the solution to the original equations will satisfy 

x + y + z — 1 = 0, 

-y + yz- z = o, 

z^ + z = 0, 

which can be easily solved to get V{I) = {(0, 1,0), (1,0,0), (1, 1, -1)}. 
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Clearly our solutions V{I) could be more complicated than isolated points. It could 
contain subsets of higher dimension. In this case, we would have a slightly different basis. 
Consider for instance the set of solutions of /i = /2 = 0. It clearly has infinitely many 
solutions, and G = {x + y + z — 1, 2y^ — 2y + 2yz + z"^ — z} is a Groebner basis. We see 
that G2 = {0} and V{l2) = C. In this case, we could not eliminate x and y, since we 
have solutions for infinitely many values of z. In summary, the number of variables we can 
actually eliminate by calculating the Groebner basis depends on the dimension of V{I). 
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